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Abstract 

We perform density functional calculations on one-dimensional zigzag 
edge graphene nano-ribbons (ZGNRs) of different widths, with and without 
edge doping including semi-local exchange-correlations. Our study reveals 
that, although the ground state of edge passivated (with hydrogen) ZGNRs 
prefers to be anti-ferromagnetic, the doping of both the edges with Boron 
atoms stabilizes the system in a ferromagnetic ground state. Both the local 
and semi-local exchange-correlations result in half-met allicity in edge pas- 
sivated ZGNRs at a finite cross-ribbon electric field. However, the ZGNR 
with Boron edges shows half-metallic behavior irrespective of the ribbon- 
width even in absence of electric field and this property sustains for any field 
strength, opening a huge possibility of applications in spintronics. 
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Nanomaterials of Carbon of different dimensionalities like Carbon nan- 
otubes, fullerenes etc. have been a subject of interest over the past few 
decades due to their potential applications in various nanoscale electronic 
devices[U El EJ. Graphene, the two dimensional flat monolayer of Carbon 
atoms packed into a honeycomb lattice is the latest addition to this family 
owing to the recent progress in experimental techniques [U El El U\- Because of 
its sophisticated low- dimensional electronic properties and huge application 
possibility, it has attracted the focus of a big scientific community in recent 
times to explore it in various aspects^ U HD1 HT 1 H21 1131 fM H5] . 

The electronic properties of nanoscale Carbon systems are governed by 
their size and geometry. Recent progress in experiments allows to make 
finite size graphene layer, termed as graphene nanoribbons (GNRs) with 
varying widths, either by cutting mechanically exfoliated graphenes[H El E] 
and patterning by electron beam lithography [16], or by controlling the epi- 
taxial growth of graphenes|17[ [T8] . Different possibilities of geometrical 
termination of graphene layer give rise to two different edge geometries, 
namely, zigzag and armchair edges, differing largely in their electronic prop- 
erties. These different edge geometries have been modeled by imposing differ- 
ent boundary conditions on Schrodinger's equation within the tight-binding 
limit [El EH1 EH E2] or on the Dirac equation for two-dimensional massless 
fermions with an effective speed of light [231 |2"H [251 [26] in previous studies. 
There exists a few many body studies also on the low-dimensional electronic 
properties [271 EH EH [30l [3T] of these systems. These also have been exten- 
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sively studied using density functional theory [3"2~| I3"3"]. 

One-dimensional antiferromagnetic zigzag edge graphene nano-ribbons 
(ZGNRs) show half-met allicity at a finite external electric field across the 
ribbon width within local density approximation (LDA) as observed by Son 
et. al. [32]. Half- metallic materials show zero band gap for electrons with one 
spin orientation and semiconducting or insulating band gap for the other, giv- 
ing rise to completely spin polarized current, which has already been observed 
in some materials like Heusler compoundsjM], manganese perovskite[35j etc. 
However, one recent study predicted that, in presence of non-local exchange- 
correlations like B3LYP, the cross-ribbon electric field makes the band gap 
of one spin orientation insulating and the other semiconducting for a fi- 
nite size edge passivated ZGNR with possible application as spin-selective 
semiconductors [36J. 

To enlighten the controversy on the half-metallicity of ZGNRs in pres- 
ence of non-local exchange-correlations, we study ID periodic ZGNRs of 
different widths with a variety of edge passivation and doping using ab ini- 
tio density functional package, SIESTA [37] . We have performed spin polar- 
ized calculations within generalized gradient approximation (GGA) consider- 
ing Perdew-Burke-Ernzerhof (PBE) exchange and correlation functionals|38] 
with double zeta polarized (DZP) basis set. GGA approximation takes into 
account the semi-local exchange-correlations which have significant impact 
on low-dimensional spin systems like GNRs. The ID ZGNRs along x-axis 
and in xy-plane are represented by the unit cell lattice vectors 2.45, 30.0, 16.0, 
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Figure 1: (Color online) Model system: a zigzag edge ID graphene with 
4 zigzag chains (4-ZGNR), translated along x-axis. The edge atoms are 
denoted by "a". The dotted arrow shows the direction of external electric 
field. 

which create sufficient vacuum, ensuring no interaction between the adjacent 
layers [32]. We have considered 400 Rydbergs energy cut-off for a real space 
mesh size and a k-point sampling of 36 k-points, uniformly positioned along 
the ID Brillouin zone. The systems considered are depicted in Fig.l with 
(z) all the edge Carbon atoms passivated by hydrogen and (ii) all the edge 
Carbons replaced (doped) with Boron (B). We vary the ribbon width (w) 
for both the cases from 4 to 10 zigzag chains, and relax all the structures for 
both ferromagnetic and antiferoomagnetic spin orientations. We then apply 
a transverse electric field along the cross-ribbon direction, i.e., along the y- 
shown in Fig.l and track the band gap variation for both the spins 
orientations. 

The undoped ZGNRs with hydrogen-passivated edges show that the an- 
tiferromagnetic ground state is stabler than the ferromagnetic state by a 
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Figure 2: (Color online) Top picture shows the accumulation of up spin and 
down spin densities on either edges of a finite segment of ID 8-ZGNR with 
hydrogen passivated edges in absence of electric field. The down panel shows 
the up spin (closed symbols) and down spin (open symbols) gaps for 6 (circle), 
8 (square) and 10-ZGNRs (triangle) as a function of electric field. Negative 
field corresponds to reverse direction. 
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small margin, of the order of meV. The ground states with zero magnetic 
moment show antiparallel alignment of spins at opposite edges as well as 
between nearest neighbors, as depicted in Fig.2. These observations are con- 
sistent with the previous LDA results [32] as well as with the theory based 
on the Hubbard Hamiltonian for bipartite lattice [39]. However, the electron- 
electron interactions can overcome such a small gap to make the ground state 
of ZGNRs ferromagnetic [27]. In absence of external electric field, both the 
upspin and downspin band gaps show insulating behavior and the insulating 
gap decreases with the increase in width, as can be seen from Fig.2. The 
gap appears to be more than the previously reported LDA gap, since LDA 
is known to underestimate the band gap. The external electric field perturbs 
the ground state band structure. Increase in field strength along positive 
direction decreases the band gap for the upspin and at a finite electric field 
strength, the up-spin channel becomes metallic. The band gap for down- 
spin however, increases with the field strength. As a result, the ID ZGNRs 
become half-metallic beyond a finite critical field strength, termed as E CjW . 
From Fig.2 it is clear that, the E C)W decreases with the increase in width (w). 
Electric field in reverse direction closes the down-spin gap and increases the 
up-spin insulating gap. Thus, different directions of external electric field 
open up different spin channels. Since, the edges have different spin po- 
larizations, the applied bias eEy (e, E and y are electronic charge, electric 
field and the distance respectively) in either directions affects the up-spin 
and down-spin band gaps in opposite fashions. Our results thus suggest that 
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in presence of gradient corrected semi-local exchange-correlations, the ID 
ZGNRs do behave as half-metallic materials beyond a certain electric field 
strength. This result strongly contradicts the observation of spin-selected 
semiconducting behavior of ZGNRs in presence of non-local exchange- cor- 
relation as found by Rudberg et al[36]. The controversy of semiconducting 
gap for one spin channel arises due to the finite size of the system, used in 
ref|36j and such a gap can be tuned with increase in length (along x-axis) 
of ZGNRs. The gradual spin density annihilation at the edges of finite size 
ZGNRs with increase in electric field as observed by Rudberg et al[3B] in 
presence of non-local exchange-correlation is nevertheless due to finite size 
effect pU] . Moreover, very high field strength as considered in ref[36] is un- 
physical (beyond perturbing regime) and leads to collapse of the structure of 
such a nano-scale system. 

The alternate idea of edge passivation is to dope the edge atoms by three 
coordinated Boron atoms, with introduction of holes in the system. We 
find that B doping at the edges makes the ground state of ZGNRs ferro- 
magnetic. Such a possibility with large hole doping in strongly correlated 
low-dimensional systems was proposed earlier [H]. In our case, the energy 
difference between the ferromagnetic and antiferromagnetic ground state is 
of the order of meV. In the ground state, the electrons at the same edges as 
well as at the opposite edges are parallely aligned (see Fig.3). In Fig.4, we 
plot the density of states (DOS) for this system for zero and nonzero elec- 
tric field. As can be seen, even in absence of any external electric field, this 
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Figure 3: (Color online) Accumulation of up spin densities on both the edges 
of a finite segment of ID 8-ZGNR with all the edge atoms doped with Boron 
in absence of electric field. 

system shows metallic behavior for the majority spin and insulating behav- 
ior for the minority spin, the typical half-metallic behavior. The increase in 
field strength in either directions does not change this behavior: the majority 
spin channel always remains conducting and the minority spin gap remains 
insulating with almost same gap for all the field strengths (see Fig. 4). It is 
because, the external field in either directions affects both the spins channels 
similarly, since the spin polarization of the opposite edges are same in this 
system. As a consequence, the change in direction of the external field cannot 
change the spin polarization of the current. The minority spin gaps for all 
the widths are much higher compared to the room temperature (0.025 eV). 
Thus, the ZGNRs with both the edges doped with B can act as a potential 
half-metallic material for any external electric field over a large temperature 
domain irrespective of its width. The projected DOS (pDOS) analysis shows 
that, the half metallicity at Fermi energy in this system is largely induced 
by the B orbitals at the edges. 

We have also studied the ZGNRs of several widths with both the edges 
doped by Nitrogen (N) atoms. This system shows antiferromagnetic ground 
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Figure 4: (Color online) Up and down spin density of states of 8-ZGNR with 
both the edges doped by Boron as a function of energy, scaled with respect 
to the Fermi energy for three different electric fields (a) — 0.2, (6)0.0 and 
(c)0.2 VA~\ 

state with different spin polarization of the opposite edges. However, al- 
though the electron accumulation on the edges is noticably low compared to 
the systems discussed above, this system shows metallic behavior for both the 
spins at Fermi energy even in absence of external electric field. The pDOS 
analysis shows that, the N orbitals contribute largely to the DOS at Fermi 
energy making the system metallic. All these observations remain consistent 
with the increase in width. 

In summary, we have studied the ID ZGNRs of different widths with edge 
passivation and with edge doping in presence of transverse electric field. Our 
study shows that, even in presence of gradient corrected semi-local exchange- 
correlations, the edge passivated ID ZGNRs can act as half-metallic materials 
beyond a certain critical field strength, which decreases with the increase in 
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ribbon width. However, Boron doped ZGNRs show half-metallic behavior 
even in absence of external electric field and this half-metallicity sustains 
for any finite field strength over a large temperature range. This observa- 
tion opens a huge possibility of application of B doped ZGNRs in spintronic 
devices. 
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